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Abstract. We calculate the diffusion coefficients of persistent random walks on cubic and 
hypercubic lattices, where the direction of a walker at a given step depends on the memory of 
one or two previous steps. These results are then applied to study a billiard model, namely a 
three-dimensional periodic Lorentz gas. The geometry of the model is studied in order to find 
the regimes in which it exhibits normal diffusion. In this regime, we calculate numerically the 
transition probabilities between cells to compare the persistent random-walk approximation 
with simulation results for the diffusion coefficient. 
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1. Introduction 

Problems dealing with the persistence of motion of tracer particles - that is, the tendency to 
continue or not in the same direction at a scattering event - are encountered in many areas 
of physics; see e.g. (H and references therein. We are specifically interested in the effect of 
persistence for the motion of random walkers on regular lattices. 

The diffusive properties of persistent random walks on two-dimensional regular lattices 
were the subject of a previous paper by two of the present authors [2 J. There, we presented 
a theory making use of the symmetries of such lattices to derive the transport coefficients of 
walks with a two-step memory. In the first part of the present paper, we extend this theory to 
hyper-cubic lattices in arbitrary dimensions, which is possible by describing the geometry of 
the lattices in a suitable way. 

Persistence effects naturally arise in the context of deterministic diffusion |f3] |4] [3] |6), 
which is concerned with the interplay between dynamical properties at the microscopic scale 
and transport properties at the macroscopic scale. A variety of different techniques are now 
available, which rely on the chaotic properties of model systems to describe their macroscopic 
properties |7), E] chap. 25]. In particular, periodic Lorentz gases and related models, such as 
multi-baker maps, are simple deterministic dynamical systems with strong chaotic properties 
which also exhibit diffusive regimes. Although the transport coefficients of these models can 
be expressed formally in terms of the microscopic dynamical properties, actually computing 
them is usually difficult, with the exception of some of the simplest toy models ll^ [T0lfTTI . 
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One reason for this is that memory effects can remain important, in spite of the chaotic 
character of the underlying dynamics. 

The diffusive properties of these models therefore provide ideal applications of the 
formalism presented in this paper. An example of this was illustrated in reference fl2l . for 
a class of two-dimensional periodic billiard tables. Extending these results, in this paper we 
apply the formalism to model the diffusive properties of higher-dimensional periodic Lorentz 
gases. 

The diffusive properties of the three-dimensional periodic Lorentz gas, which consists 
of the free motion of independent tracer particles in a cubic array of spherical obstacles, are 
interesting in their own right. In two spatial dimensions, the existence of diffusive regimes 
in such systems has been rigorously established lfT3l H4l . It relies on the finite-horizon 
property, which requires that the system admits no ballistic trajectories, i.e. those which never 
collide with any obstacle. In this case, it is possible to change scales from microscopic to 
macroscopic, reducing the complicated motion of tracer particles at the microscopic level to 
a diffusive equation at the macroscopic level. When the horizon is infinite on the other hand, 
there is rather a weakly superdiffusive process, with mean-squared displacement growing like 
t log? 03] [16), as recently shown rigorously in 1171 . 

The necessity of finite horizon to have normal diffusion in two dimensions led to the idea 
that this was also necessary in three dimensions - see, for example, reference 1 18 1. Recently, 
however, it was argued by one of the present authors 1 19 1 that in higher-dimensional billiards, 
normal diffusion, by which we mean an asymptotically linear growth in time of the mean- 
squared displacement, may arise even in the absence of finite horizon. In fact, three different 
types of horizon can be identified in the three-dimensional periodic Lorentz gas. The key 
observation is that it is only "planar" gaps - those with infinite extension in two dimensions 
- which induce anomalous diffusion. If there are only "cylindrical" gaps, whose extension is 
limited to a single dimension, then the available space in which particles can move ballistically 
is limited. This leads to a decay of correlations which is fast enough to give normal diffusion 
at the level of the mean-squared displacement, although higher moments of the displacement 
distribution may be non-Gaussian [i T9ll . 

The paper is organized as follows. Section [2] describes the computation of the transport 
coefficient of walks on hypercubic lattices with one and two-step memories. In the second 
part of this paper, we apply this formalism to the diffusive regimes of the three-dimensional 
periodic Lorentz gas. In section [3] we give a detailed description of the three-dimensional 
periodic Lorentz gas introduced in reference [19], in particular delimiting the regimes with 
qualitatively different behaviour in parameter space. We then apply the results on persistent 
random walks to the diffusive regimes of this model in section [4] Conclusions are drawn in 
section |5] 

2. Persistent random walks on cubic lattices 

In this section, we describe a way to incorporate the specific geometry of cubic and hyper- 
cubic lattices in the framework presented in reference 1 2 1 for calculating diffusion coefficients 
for persistent random walks on lattices. 

We start by considering the motion of independent walkers on a regular cubic lattice in 
three dimensions. Given their initial position vq at time t = 0, the walkers' trajectories are 
specified by the sequence {vq, ■ ■ ■ ,v„} of their successive displacements. Here we consider 
dynamics in discrete time, so that the time sequences are simply assumed to be incremented by 
identical time steps T as the walkers move from site to site. In the sequel we will loosely refer 
to the displacement vectors as velocity vectors; they are in fact dimensionless unit vectors. 
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The sequence of successive displacements is determined by the underlying dynamics, 
whether deterministic or stochastic. At the coarse level of description of the lattice dynamics, 
this is interpreted as a persistent type of random walk, where some memory effects are 
accounted for: the probability that the nth step is taken in the direction v„ depends on the 
past history v n -\,v n -2, ■ ■ ■■ 

The quantity of interest here is the diffusion coefficient D of such persistent processes, 
which measures the linear growth in time of the mean-squared displacement of walkers. This 
can be written in terms of velocity autocorrelations using the Taylor-Green-Kubo expression: 



D 



2dr 



1 



2 lim 



Ic 

E 

«=i 



VQ ■ V„ 



(2.1) 



where d denotes the dimensionality of the lattice, here d = 3, and I is the lattice spacing. The 
(dimensionless) velocity autocorrelations are computed as averages (•) over the equilibrium 
distribution, denoted by ju, of the underlying process, so that the problem reduces to 
computing the quantities 

(v -v n )= £ v -v n n({v Ql . ..,«„}). (2.2) 

«« 

Following the approach of reference Q, we wish to compute the terms in this sum, and 
hence the corresponding diffusion coefficient ( |2.1[ ), for three different types of random walks, 
namely those with zero-step, single-step and two-step memories. These cases all involve 
factorisations of the measure fi({vo, . . ■ , v„}) into products of probability measures which 
depend on a number of velocity vectors, equal to the number of steps of memory of the 
walkers. These measures will be denoted by p throughout. 

The schemes we outline below allow to write equation (2.2 1 as a sum of powers of 
matrices, so that ( 2. 1 1 boils down to a geometric series, which can then be resummed to obtain 
an expression for the diffusion coefficient that is readily computable given the probabilities 
that characterise the allowed transitions in the process. 



2.1. Description of geometry of cubic lattices 

It is first necessary to find a succinct description of the geometry of the cubic lattices that 
we wish to study. The six directions of the three-dimensional cubic lattice and corresponding 
displacement vectors are specified in terms of the unit vectors e, of a Cartesian coordinate 
system as ±e,-, i— 1,2,3. 

The crucial property required for the application of our method is that all of these unit 
vectors can be obtained by repeated application of a single transformation G, which generates 
the cyclic group 

& = {G; = G', / = 0,...,5}. (2.3) 
One possible choice of G gives the following group elements: 

0-1 

1 
v 1 

G 5 = G 2 = [ -II. (2.4) 



-G = G 3 




Diffusive properties of persistent walks on cubic lattices 



4 



Figure [T] displays the six possible directions of a walker on this lattice, numbered 
according to repeated iterations by G. Thus a walker with incoming direction e\, indicated 
by the arrow, can be deflected to any of the six directions G'ei, i = 0,...,5, corresponding 
respectively to ei, e2, 63, — ei, — e?, and — ej,. 



2 




Figure 1. The possible directions of motion on a cubic lattice, labelled from to 5 relative 
to the incoming direction shown by the arrow. These directions are obtained by successive 
applications of the transformation G given in equation (2.4). 



A similar transformation G can easily be identified for a walk on a of-dimensional hyper- 
cubic lattice: 



(2.5) 



which maps the unit vectors onto the 2ii-cycle e\ n- eo i-> • • • n- e f / n- — e\ i->- • • • i->- — e ( (. 
2.2. No-Memory Approximation (NMA) 
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We now proceed to calculate the diffusion coefficient (2.1 1 for random walks with different 
memory lengths. The simplest case is that of a Bernoulli process for the velocity trials, so 
that the walkers have no memory of their history as they proceed to their next position. The 
probability measure jj, thus factorises: 



H{{v ,...,v„}) =Y[p(vi). 



(2.6) 



Given that the lattice is rotation invariant and that p is uniform, the velocity autocorrelation 
(12. 2b vanishes: 



{vo-v n ) 



(2.7) 
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The diffusion coefficient of the random walk without memory is then given by 

e 2 

Dnma — ^rr ■ 
2ax 



(2.8) 



2.3. One-Step Memory Approximation (1-SMA) 

We now assume that the velocity vectors obey a Markov process, for which v„ takes on 
different values according to the velocity at the previous step v n ~\. We may then write 



H({vq, v„}) = n^Mwi-iMwo) ■ 



(2.9) 



Here, P(v'\v) denotes the one-step conditional probability that the walker moves with 
displacement v', given that it made a displacement v at the previous step. 

Considering for definiteness the three-dimensional lattice and using the elements of the 
group Sf, we express each velocity vector Vk in terms of the first one, vq, as Vk = G' k vo, where 
each if. g {0, ...,5}. Substituting this into the expression for the velocity autocorrelation 



vq ■ v„), equation ( |2.2| , we obtain, using the factorisation (2.9 1, 

i)/>Oo) = 



v -v„Y\P(vi\vi 

VQ,...,V n i=\ 

In this expression, 



6 

£ 



V o • G M V m;; , 



m 



;f > (G'"t;o|G , "- 1 t;o) 



(2.10) 



(2.11) 



are the elements of the stochastic matrix M of the Markov chain associated to the persistent 
random walk, and p,- = p(e,-) are the elements of its invariant (equilibrium) distribution, 
denoted P, evaluated with a velocity in the z'th lattice direction. The invariance of P is 
expressed as m i,jPj — Pi- The same notations were used in Q and will be used throughout 
this article. 

The terms involving M in (2.10 1 constitute the matrix product of n copies of M. 
Furthermore, since the invariant distribution is uniform over the lattice directions, we can 
choose an arbitrary direction for vq, and hence write 



: v -v m\ [+v 



r (") 



+ ■ 



r-5 W 



(2.12) 



where m-y denote the elements of M". 

The actual value of the diffusion coefficient depends on the probabilities P(G^v\v), 
which are parameters of the model, subject to the constraints Y l jP(G J v\v) = 1. To simplify 
the notation, we assume rotational invariance of the process, i.e. independence with respect to 
the value of v, and we denote the conditional probabilities of these walks by Pj = P(G^v\v), 
where j = 0, . . .,5. 

The transition matrix M given by ( |2.11| > is thus the cyclic matrix 

/ P Pi Pi Pi Pa Ps \ 

Ps Pq Pi Pi Pi A 

Pa P5 Po Pi Pi A 
M = (2 13) 

Pi Pa P 5 Pq A ft 1 v 

ft ft ft ft ft ft 

V ft ft ft ft ft ft / 
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The matrix M" shares the same property of cyclicity, so that it also has only six distinct entries. 
It is thus possible to proceed along the lines described in [2 1 and obtain the recurrence relation 



/ 


In) 
m U 


— m 




(«) 
m 2,l 


— m 


V 


(«) 
m 3,l 


-m 



P) 


-Pi Pi 


-Pa Pi 


-P5 


P5 


-Pi Pq 


-P3 Pi 


-Pa 


Pa 


-Pi P5 


-Pi Pq 


-Pi 


Pq 


-P 3 Pi 


-Pa Pi 


-P5 


P5 


-Pi Pq 


-Pi Pi 


-Pa 


P4 


-Pi P5 


-Pi Pq 


-Pi 



(n-l) (h-1) 
m l,l " m 4,l 




(2.14) 



[Note that the left-hand side of this equation was chosen to reduce the size of the matrix 
involved and to calculate the element required in ( |2.12| >.] As a consequence, we can write for 
the velocity autocorrelation (2.12 1 



P0-P3 Pi -Pa P1-P5 \ " 1 / ^0-^3 
va v„) (1 0) I P5 -Pi Pq- Pi Pi - Pa Pi-Pa 
Pa-Pi P5-P1 P Q -Pi J \ P2-P5 

and thus obtain the expression of the diffusion coefficient ( |2.1| i as 

-1 



(2.15) 



^lSMA 
^NMA 



1+2 1 



i+Pi-Ph P4-P1 P5-P2 
P2-P5 1+P3-P0 P4-P1 
Pi-Pa P2-P5 y+P}-Pb 




(2.16) 



by using the result that £~ = q A" = ( I — A) ~ 1 , where I is the identity matrix, for a square matrix 
A whose eigenvalues are all strictly less than 1 in modulus, 

This result easily generalises to a hyper-cubic lattice in any dimension d. Note also that 
for a symmetric process, in which Pi = P4 and P% = P5, we recover the diffusion coefficient 

I+B0-P3 
l-Po+Pi' 
in agreement with the result stated in 0. 



D 



ISM A 



f NMA : 



(2.17) 



2.4. Two-Step Memory Approximation (2-SMA) 

Let us now suppose that the velocity vectors obey a random process for which the probability 
of v n takes on different values according to the velocities at the two previous steps, t)„_ 1 and 
v n -2, so that we may write 

n 

H{{v , ■ ■ .,«„}) = YY p { v i\ v i-ii v i-l)p{ v Q, v l) ■ ( 2 - 18 ) 

1=2 

The velocity autocorrelation \2.2\ function is then 

(2.19) 

Since the probability transitions P(vj\vi- i,Uj_2) have symmetries similar to those used 



(vo-v n )= Y, v -v„Y[P{vi\vi-i,Vi^ 2 )p(vQ,vi) . 

{v, v } 1=2 



in reference [2 |, the computation of equation (2.19 1 reduces to an expression very similar to 
that found there for walks on one- and two-dimensional lattices. The details of the derivation 
are a bit more involved than the one-step memory persistent walks, so we will limit ourselves 
to stating the results. 
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Letting z = 2d denote the coordination number of the lattice, and writing] Pj^ = 
P(G z ~ k G z ~-iv\G z ~jv 7 v), which is the conditional probability of making a displacement v 
given that the two preceding displacements were successively G-'G v and G k v, we define the 
Z x z matrix 



K(<j>) 



I 



\ 



PqO 
<t>Poi 



Ml 



^-1,0 



(2.20) 



The argument (j) in this expression is a complex number such that (j) z = 1. In the case of 
two-dimensional lattices, only two of these roots are relevant, corresponding to the complex 
exponential of the smallest angle between two lattice vectors, = exp(±2/7r/z). For hyper- 
cubic lattices in arbitrary dimensions, however, we must consider a priori all the z possible 
roots of unity, (j)j = e\p(2inj /z), j = 0,...,z— 1. 

A direct calculation of ( |2.19 1 shows that the velocity autocorrelation takes the form 



(vo ■ v„) 



l-.-l 



--1 



£ ai ^r 1 diag(l,0 / ,...,0f 



Pi 



(2.21) 



where diag ( 1 , ; , . . . , 0? 1 ) denotes the matrix with elements listed on the main diagonal and 
elsewhere. For the three-dimensional cubic lattice, the coefficients aj are found to be 



fl() = fl2 = a A = 0, 

a\ = «3 = «5 = 2, 



(2.22) 



which compares to a\ = «3 = 2 and oq = a^ = in the case of the two-dimensional square 
lattice [2]. In the case of a af-dimensional hyper-cubic lattice, this generalises to 



(2.23) 



a 2j = 0, j = 0, . . . ,d— 1, 
02.7+1=2, j = 0,...,d-l, 
The diffusion coefficient of a two-step memory persistent random walk on a d- 



dimensional hyper-cubic lattice is thus 



O2SMA 



£*NMA 



= 1+4 l-.-l <^ "EVz-Kifoj-iT'&agihfr 



1j-U 



J =1 



&J-1J 



Pi 



Pz 



(2.24) 



where /, denotes the z x z identity matrix. 



3. Three-dimensional periodic Lorentz gas 



Equations (2.8 1, (2.16 1 and (2.24 1 can be put to the test to probe the diffusive regimes of 
periodic Lorentz gases. The diffusive motion of the tracers results from the chaotic nature 
of the microscopic dynamics and the fast decay of correlations, which are in turn due to 
the convex nature of the obstacles. Taking into consideration the different diffusive regimes 
of these models, which, as we argued earlier, depend on the nature of their horizon, we 



% This expression differs from that given in [ 2 1 due to a typographical error in that paper - they are really the same. 
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investigate how the microscopic dynamical properties of the system determine the diffusion 
coefficient. 

Machta and Zwanzig [20] addressed this issue in a particular limiting case, showing that, 
in the limit where the obstacles are so close together that a tracer will remain localised on 
ekach lattice site for a very long time (compared to the mean time separating two collision 
events), the process of diffusion on the Lorentz gas is well approximated by the dimensional 
prediction (2.8 1, where the lattice spacing £ is the distance separating two neighbouring 
obstacles and z is the trapping time, which can be computed in terms of the geometry of 
the billiard as a simple consequence of ergodicity. That is to say, when the geometry of the 
billiard is such that two neighbouring disks nearly touch, the Lorentz gas is well approximated 
by a Bernoulli process, modeling the random hopping of tracers from cell to cell, with time- 
and length-scales specified according to the geometry of the billiard. 

Different approximation schemes have been proposed to go beyond this zeroth-order 
approximation and account for corrections to it GT1 l22l ; see, in particular, reference 1231 for 
an overview. A consistent approach to understanding the effect of these corrections in two- 
dimensional diffusive billiards was described in lfT2l . The idea is to approximate the hopping 
process of tracer particles by persistent random walks with finite memory, and thus estimate 
the diffusion coefficient of the billiard by the two-dimensional lattice equivalents of the one- 
or two-step formulas ( 2. 16 1 and ( 2.24 1. 

We discuss below the transposition of these results to the diffusive regimes of the three- 
dimensional periodic Lorentz gas. 



3.1. Geometry of simple three-dimensional periodic Lorentz gas model 

We begin with a detailed description of the geometry and the different horizon regimes of the 
system studied in reference fl9l ; additional details are given in reference ll24ll . 

The model consists of a three-dimensional (3D) periodic Lorentz gas constructed out of 
cubic unit cells of side length I, having eight "outer" spheres of radius p out £ at its corners 
and a single "inner" sphere of radius p m l at its centre - see figure [2] The infinitely-extended 
periodic structure formed in this way is symmetric under interchange of p m and p out ; without 
loss of generality, we take p ou t > Pin- 




Figure 2. Geometry of the obstacles in a single cell of the 3D periodic Lorentz gas model for 
Pout — 0.45£ and Pj n = 0.3(M, in the cylindrical-horizon regime. 

This model seems to be the simplest one which allows a finite horizon, although this 
is possible only when the spheres are permitted to overlap. It is known that finite-horizon 
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periodic Lorentz gases with non-overlapping spheres in fact exist in any dimension |25|, but 
we are not aware of any explicit constructions of such models, even in the case of three 
dimensions. 



Lattice of outer spheres The spheres of radius p out £ form a simple cubic lattice. This lattice 
has the following properties: 

• When p ou t < 1 /2, the spheres are disjoint. In this case, there are free planes &5l 
in the structure, that is, infinite planes which do not intersect any of the spheres, in 
particular there are free planes centered on the faces of the unit cell. In this case, we 
say that there is a planar horizon (PH). When p ou t is small, there are additional planes at 
different diagonal angles, analogously to the two-dimensional infinite-horizon Lorentz 

gas oa nana. 

• When p ou t > 1/2, the spheres overlap, thereby automatically blocking all planes. The 
overlaps (intersections) of the spheres partially cover the faces of the cubes, leaving a 
space in between which acts as an exit towards the adjacent cell. 

• When Pout > 1 /V2> the overlaps completely cover the faces of the unit cell, so that it is 
no longer possible to exit the cell. 

• When Po U t > v3/2, all of space is covered, and it is no longer possible to define a billiard 
dynamics. 



Conditions for normal diffusion: cylindrical horizon As shown in reference [19|, the 
necessary and sufficient condition to have normal diffusion is that all free planes are blocked; 
if there are free planes, then the diffusion is weakly anomalous. The conditions to block all 
planes are as follows. 

• All free planes are automatically blocked for p ou t > 1 /2, when the p ou t-spheres overlap. 

• If the Pout-spheres do not overlap, then it is necessary to introduce the p; n -sphere to block 
planes which are parallel to the faces of the unit cell. For this blocking to occur, we need 

Pin > 1/2 - Pout- 

• Furthermore, we must also block diagonal planes at 45 degree angles, which requires 
that p om > 1/(2^/2) or p in > 1/(2^2). 

If all of these conditions are satisfied, then we no longer have free planes, but may have free 
cylinders ("cylindrical gaps") in the structure; we then say that there is a cylindrical horizon 
(CH). 



Conditions for finite horizon Stronger statistical properties - e.g. faster decay of correlations 
- may be expected when there is a finite horizon 11811191 , i.e. where the length of free paths 
between collisions with obstacles is bounded above. To obtain this, not only all planar gaps, 
but also all cylindrical gaps must be blocked, i.e. all holes viewed from any direction must be 
blocked. To do so, the following conditions must be fulfilled: 

• The Pout-spheres must overlap, p ou t > 1/2. Furthermore, the projection of the p; n - 
sphere on each face of the unit cell must cover the available exit space, as illustrated 



in figure 3(a) Letting d be the maximum width of overlap of the resulting discs of radius 
Pout on a face of the unit cell, we have d 2 = p 2 ut — 1/4, and we need p; n > 1/2 — d to 
block the space. 
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Pout / /\ 

/ 1 I 




(a) Face of unit cell 



(b) Mid-plane of unit cell 



Figure 3. Geometry of the 3D periodic Lorentz gas. (a) Cross-section of the unit cell in one 
of its faces. The overlapping outer spheres of radius p out > 1/2, give rise to four overlapping 
discs (shown in green); the maximum width of their overlap is denoted d. The central disc 
(red) shows the minimum radius p- m = 1/2 — d of the central sphere such that its projection 
covers the gap between the p out -discs on the face, (b) Geometry of the mid-plane of a unit cell 
for parameters giving a finite horizon. The outer discs are cross-sections of the overlaps of the 
outer p 011 t-spheres, and have radius d equal to the overlap parameter in (a). The inner disc is 
the cross-section of the inner Pj n -sphere. 




Figure 4. Finite horizon can be achieved in a three-dimensional lattice, provided the spheres 
are allowed to overlap. Here the available space for diffusing particles is shown for parameter 
values p olll = 0.651 and p m = Q.\5( (in the FH1 region) in an unfolded channel. 



• We must block cylindrical corridors which cross the structure at a 45 degree angle at 
the level of the mid-plane of a unit cell, which corresponds to the planar cross-section 
with most available space in the unit cell. The mid-plane has the geometry shown in 
figure 3(b) with four outer discs of radius d, and a central disc of radius p m ; these discs 
are the intersection of the p out -overlaps and of the p; n -sphere, respectively, with the mid- 
plane. Free diagonal trajectories in this plane at an angle of 45 degrees give rise to small 
cylindrical corridors. These will be blocked if there is no free line in the mid-plane. This 
blocking occurs provided either d>l/ (2\/2), i.e. p out > V3/(2V2, or if p in > \/(2y/2), 
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thus giving rise to two distinct finite horizon regimes (FH1 and FH2), which are in fact 
disjoint. 

Figure [4] depicts the space available for tracer particles in a channel of three consecutive cells 
for a particular finite-horizon case. 

Localisation of trajectories Having fixed p out , it is also necessary to calculate the value of 
Pi n above which the trajectories become localised (L) between neighbouring spheres, and 
are thus no longer able to diffuse. For p out < l/v2, when there are still exits available 
on the faces of the cubic unit cell, this happens exactly when the discs in the mid-plane 
touch, i.e. when p nl +d = 1/ v2, so that the condition for localised trajectories becomes lfl9ll 
Pin> l/%/2-v/pi ut -l/4. 



Condition to fill space Finally, we calculate when the spheres fill all space (denoted U, for 
undefined): 

• When p ou t < l/v2, this occurs when the pj n -spheres are large enough that their 
intersection with each face of the cube, which is a disc of radius yj p? — 1 /4, covers 
the exit on a face left open by the p out -spheres. This gives the condition p? n > p£ ut + 
l/4-VPo 2 ut-l/4. 

• When p ou t > 1 / \/2, the condition is that p; n be large enough to cover the space left by 
the Pout-spheres inside the unit cell. The condition can again be found by looking at 
the mid-plane, where there is most available space: the disc of radius p; n must cover 
the space left by the discs of radius d (which are cross-sections of the overlaps of the 
Pout-spheres). This occurs when p; n > 1/2 — \J p^ ut — 1 /2. 



Parameter space The complete parameter space of this model is shown in figure [5] 
exhibiting the regions in parameter space corresponding to the regimes of qualitatively 
different behaviour discussed abov^§] Note that if pj n > V3/2 — p ont , then the p out - and p m - 
spheres overlap, and if p m > 0.5 then neighbouring p; n -spheres also overlap. These conditions 
are marked by the dotted lines in the figure. 



4. Persistence in the diffusive regimes of the three-dimensional Lorentz gas 

In this section, we study the dependence of the diffusion coefficient on the geometrical 
parameters of the 3D periodic Lorentz gas model in the finite- (FH1) and cylindrical-horizon 
(CH) regimes, comparing the numerical results with the finite-memory approximations (2.8 i, 
( |2~T6l ) and $22ty . 



4.1. Approximation by the NMA process 

The computation of the dimensional formula ( |2.8| l relies on that of the residence time T. An 
exact formula is available for this quantity |27T 

. \Q\ \s 2 \ 



\dQ\ |BT 



(4.1) 



§ A similar diagram of parameter space for a two-dimensional version of the model was given in reference [ 26]. 
However, the symmetry between p ollt and p m was overlooked there; see also reference [19]. 
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Figure 5. Parameter space of the three-dimensional periodic Lorentz gas as a function of 
the geometrical parameters p ou t and p m . Solid lines divide regimes of qualitatively different 
behaviour, which are also shaded with different colours and labelled as follows: PH: planar 
horizon; CH: cylindrical horizon; FH1 and FH2: finite horizon; L: localised, non-diffusive 
motion; U: undefined (all space filled). Note that the FH regime is divided into two disjoint 
regions. The dashed lines mark the different conditions referred to in the text. The diagonal 
dotted line separates regions where the pj n -spheres do (above) and do not (below) overlap the 
p oll t-spheres. The diagram is reflection-symmetric in the line p; n = p out , but for clarity only 
the lower half is shown. 



where \Q\ denotes the volume of the billiard domain outside the obstacles, \dQ\ the surface 
area of the available gaps separating neighbouring cells, \S 2 \ — 4-n the surface area of the unit 
sphere in three dimensions, and \B 2 \ =n the volume (area) of the unit disk in two dimensions, 
and we assume unit velocity. The explicit formulas giving the values of \Q\ and \ dQ\ are rather 
lengthy and will not be given here; see reference [28 1. 

The validity of equation ( |4. 1[ ) can be tested by comparison with numerical computation 
of the residence time, as shown in figure [6] Here, and in the remainder of the paper, we 
restrict attention to values of p out close to the limiting value 1 / v2 and p m close to 0, so that 
the geometry is that of a single, cubic unit cell. 

4.2. Approximation by the 1SMA and 2SMA processes 

Single- and two-step memory processes can be derived as approximations, at the lattice level, 
to the dynamics of the Lorentz gas. This is done by computing numerically the statistics 
of tracer particles as they jump from cell to cell, so as to estimate the single- and two-step 
memory probability transitions. 

The results are shown in figure [7] for the single-step memory process, where the six 
transition probabilities Pj, i = 0,...,5, are displayed as functions of the outer radius p out for 
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0.00 0.05 0.10 0.15 0.20 
6 = l/V2"-Pout 



Figure 6. Residence time X, equation \4. l\ , compared to direct numerical simulations. The 
results are shown for two values of the inner radius, p m = and p m = 0.15, as functions of 
5 = \j\fl — p out , which is the characteristic size of the gaps separating neighbouring cells. 
The curves are very similar since the volume of the inner sphere remains small. In this and the 
following results we take t. = \. 



different values of the inner radius p m . 

For the two-step process, the computation of the transition probabilities Pi j is shown in 
figure[8]for p m = 0, that is in the absence of a sphere at the center of the cell. The six different 
panels each correspond to a given i = 0, . . . ,5. The same is shown in figure|9]for pj n =0.15. 





0.00 0.05 0.10 0.15 0.20 0.00 0.05 0.10 0.15 0.20 0.00 0.05 0.10 0.15 0.20 

<*> = l/VT-p out ii=l/VT-p out rf= l/V2"-p ou , 

(d) (e) (f) 



Figure 7. Numerical computations of the probabilities Pq , . . . , P5 of the single step memory 
process, appearing in l |2.13) . The six panels shown correspond to as many different values of 
Pin, where the probabilities are shown as functions of 5. The dashed line at = 1/6 indicates 
the value for a memory less (NMA) walk. Here and in figures [8] and [9] the conventions are 
as follows: Empty squares (blue), Pq; empty upward triangles (cyan), Pi; empty downward 
triangles (green), P2; filled squares (red) P3; filled upward triangles (magenta), P4; filled 
downward triangles (brown), P5. In all cases we verify the symmetry Pi = Pj = P4 = P5, 
which also remain close to 1/6. 
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Figure 8. Numerical computations of the 36 probabilities P, j which appear in |2.20| , 
corresponding to a cell with no sphere at its center, i.e. the inner radius p m = 0. The symmetries 
of the process are reflected by the similarities between figures [8(b)] |8(c)||8(e)| and |8(f)| The 
colour coding is similar to figure|7] 




Figure 9. Numerical computations of the probabilities 36 Py which appear in |2.20| , 
corresponding to inner radius p m = 0.15. Here again the symmetries of the process are 
reflected by the similarities between figures |9(b)||9(c)||9(e)| and |9(f)| 
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4.3. Diffusion coefficient of the billiard 




Figure 1 0. Di ffusion coefficient normalised with respect to the dimensional prediction Z>nma, 
equation 



2.S 



, vs. gap size S = l/v2 — p oul , plotted for different values of the inner radius 
Pin = 0, ...,0.15. The symbols (black) correspond to direct numerical computation of this 
quantity, the long dashed (green) lines to the single-step memory diffusion coefficient j2.16( , 
and the solid (red) lines to the two-step memory diffusion coefficient \2.24) . The vertical 
dotted lines indicate the separation between the finite and infinite horizon regimes. 



Having computed the probability transitions associated to the single and two-step 
memory processes, we can compute the invariant distribution P and substitute the results into 
equations (2.16 1 and ( 2.24 1 to obtain values of the diffusion coefficients. These are compared 
to the diffusion coefficient of the billiard calculated from direct simulations in figure [TO] 

We can draw several conclusions from the results shown in figure 10 Firstly, we remark 
that in the 3D model studied here there is relatively little back-scattering, i.e. motion in 
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which the particle reverses its direction between arriving and leaving a given cell. This gives 
an important contribution to the diffusion coefficient, and, in particular, corresponds to the 
fact that here we find that the diffusion coefficient is larger than the memoryless (NMA) 
approximation, while in reference [12] the diffusion coefficient tended to lie below the results 
of this approximation. Note, however, that this effect depends strongly on the particular model 
used. 

It is also interesting to note that in the finite-horizon regime, i.e. left of the dotted vertical 
lines in figures 10(b)p0(f) approximating the diffusion coefficient by the one-step memory 
process (2.16 1 is just as good as the two-step process ( 2.24 1. In the cylindrical-horizon regime, 
however, the two results are different; the single-step approximation gets poorer as p ou t 
decreases, whereas the two-step process yields more accurate estimates. This corresponds 
to the fact that correlations decay more slowly in the cylindrical-horizon regime lfT9l . so that 
memory effects persist for longer. 



5. Conclusions 



The cyclic structures of certain regular lattices underly symmetries of their statistical 
properties which can be exploited to greatly simplify their analysis. Examples are two- 
dimensional lattices such as the square, the honeycomb and the triangular lattice, which 
were studied in reference [2]. Other examples include, in higher dimensions, the hypercubic 
lattices studied in this paper. Having exhibited the cyclic structures of these lattices, we were 
able to extend our previous results to hypercubic lattices with suitable adaptations, in order to 
calculate the diffusion coefficients of persistent random walks with up to two steps of memory. 

Our method is especially useful to compute the correlations of persistent walks on such 
regular lattices. In particular, the velocity autocorrelations of a two-step persistent walk may 
be recast in terms of matrix powers, which can then easily be resummed to obtain a readily- 
computable expression for the diffusion coefficient. 

Among the many applications of persistent random walks, deterministic diffusive 
processes are ideal candidates to apply our method. The three-dimensional periodic Lorentz 
gas is particularly interesting as it exhibits two distinct types of diffusive regimes, one with 
finite horizon, where memory effects decay fast, and another with cylindrical horizon, where 
memory effects can remain important. In this latter case, the approximation of the diffusive 
process by a two-step memory walk proves much more accurate than the single-step process. 

We remark that the application of our formalism to the diffusive properties of Lorentz 
gases relies on the numerical computation of the transition probabilities corresponding to 
the persistent process with which we approximate the deterministic process. Since there 
are 30 transition probabilities for the two-step memory walk, their analytical calculation is 
a daunting task. It relies on knowledge of the statistics of trapped trajectories and involves 
contributions from different time scales. Nonetheless, this computation is formally possible, 
and is in principle much simpler than that of the actual diffusion coefficient. 
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